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Abstract. A quantitative prediction of Conformal Field Theory (CFT), which relates 
the second moment of the energy-density correlator away from criticality to the value 
of the central charge, is verified in the sine-Gordon model. By exploiting the boson- 
fermion duality of two-dimensional field theories, this result also allows to show the 
validity of the prediction in the strong coupling regime of the Thirring model. 



PACS numbers: 11.10, 11.25.H 



§ To whom correspondence should be addressed (naon@fisica.unlp.edu.ar) 



On a CFT prediction 



2 



Some time ago Cardy U derived a quantitative prediction of conformal invariance |2|. [J 
for 2D systems in the scaling regime, away from the critical point. Starting from the so 
called 'c-theorem' [|J, he was able to relate the value of the conformal anomaly c, which 
characterizes the model at the critical point, to the second moment of the energy-density 
correlator in the non-critical theory: 

A tfW«' iri! (2-A,)" (1) 
where e is the energy- density operator, A E is its scaling dimension and t oc (T — T c ) is the 
coupling constant of the interaction term that takes the system away from criticality. It 
is interesting to notice that a similar sum rule has been recently obtained by Jancovici || 
in the context of Classical Statistical Mechanics. This author considered the correlations 
of the number-density of particles in a 2D two-component plasma (Coulomb gas). 

The validity of ([!]) has been explicitly verified for the Ising model 0, for 2D self- 
avoiding rings ||, and for the Baxter model 0. In this last case the formula could be 
checked only in the weak-coupling limit, by describing the system in terms of a massive 
Thirring model jS| and performing a first-order perturbative computation. 

The main purpose of this comment is to show that ([I]) also holds for a bosonic 
QFT with highly non-trivial interactions, the well-known sine-Gordon (SG) model with 
Euclidean Lagrangian density given by 

L = i(9 M $) 2 -^cos(v / A«I») + p (2) 

where a and A are real constants. 

In this work we shall perform a perturbative computation up to second order in A. 
We are then naturally led to consider the renormalization of this theory. Fortunately 
this issue has been already analyzed by many authors || [R], [11], 12, 13, 14]. One of 
the main conclusions is that for A < 8tt a normal order procedure that eliminates the 
contributions of the tadpoles is enough to have a finite theory. The only effect of this 
prescription is to renormalize the constant a. We shall then restrict our study to this 
case. 

Since we want to verify equation (|I|) for the model given by (|2|), we will take 
e = (1 — cosV\&)/\ and t = a. Taking into account that the model of free massless 
scalars has a conformal charge c = 1 and that the scaling dimension of e for this case is 
equal to A/47T, (H) reads 

F(a, A) = / d 2 x \x\ 2 (e(x)e(O)) = -± 3—. (3) 
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Expanding the interaction term up to order A 2 one has e = |$ 2 — + ^r$ e 
Replacing this expression in (^) we obtain 

F(a, A) = A(a, A) + B(a, A) + C(a, A) + D(a, A), 

where 

A(a, *) = \f d 2 x \x\ 2 ($ 2 (x)$ 2 (0)) Q , 
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B(a, A) = - — J d 2 x \x\ 2 ($ 2 (x)$ 4 (0)) Q , 
C(a, A) = J d 2 x \x\ 2 ($ 4 (0)$ 4 (x)) Q 



D{a, A) = ^ / d 2 x |x| 2 ($ 2 (x)$ s (0)) Q . 



and 

6! 

At this point we notice that D(a, X), up to this order, contains only tadpoles which, 
as explained above, were already considered in the renormalization of a. Then we must 
disregard this contribution in the present context. Now, in order to illustrate the main 
features of the computation, we shall briefly describe the evaluation of A(a,X), which 
involves both analytical and numerical procedures. In the above equations ( ) a means 
v.e.v. with respect to the SG Lagrangian expanded up to second order in A. From now 
on we will decompose this Lagrangian into free and interaction pieces as 

^o = ^$) 2 + f$ 2 , 



aX ^ A aX 2 

J int — 



$ 4 4- $ fc 



4! 6! 

Using Wick's theorem and the well-known expression for the free bosonic propagator 
($(x)$(0)) = l/(2n)K Q (- s /a\x\) (K is a modified Bessel function of zeroth order), after 
a convenient rescaling of the form x — > -7=, we obtain 

A{a,X) = 2a2 ^ 2n y I d 2 x\x\ 2 Kl{\x\) + ^ X J J d 2 xd 2 Xl \x\ 2 K 2 (\ Xl - x\)K 2 (\ Xl \) 



X 2 

+ 



a 2 (2n) 



1 1 J d 2x d 2 xi d 2 x 2 \x\ 2 \^Kl{\x 2 - x 1 \)Kq(\x 1 \)Kq(\x 2 - x\) 



+ \kI{\x 2 - Xl \)K Q [\x 2 \)K [\x x - x\)K (\x\) 
o 

+ -K$(\x 2 -xi|)if (M)-M|z2 - xDKoi^DKadxt - x\)}. 

The first two terms are related to tabulated integrals [15| yielding the first order 

analytical expression F(a, A) = (1 + A/(47t))/(127tq; 2 ) which can be shown to verify (|3|) 

in a straightforward way. Concerning the remaining second order terms we have three 

contributions corresponding to the prefactors 1/8, 1/6 and 1/4 in the last integrand. 

Due to its symmetry, the first one can be analytically computed by repeatedly using 

r 2n 
J d 2 x\x\ 2 K 2 (\x-y\) = — + n\y\ 2 . 

The result is 2ir 3 . The computation of the other two contributions is more involved. 
In fact we could not find analytical results in these cases. However we were able 
to considerably simplify these multiple integrals in order to facilitate their numerical 
evaluation. Indeed, using the Fourier transform of Kq and the integral representation 
of the first kind Bessel function J : 

Jq(\x\) — — I d9 exp (— i\x\ cos6) , 

Z7T JO 
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xK (\x\) = 4:(2n) 3 J dr J dp 



we obtain 

d 2 xd 2 xi d 2 x 2 \x\ 2 Kq(\x2 — Xi\)K (\xi\)K (\x 2 — x\) 

rp (1 — p 2 )Kq (r) Jo{pr) 
(p 2 + l) 5 ' 

Using NIntegrate in the program Mathematica for the double integral in this expression 
one obtains the value 0, 04874. Similar manipulations with the last contribution to 
A(a, A) led us to the following numerical results: 

drjdp J dk {p2 - i)4(fc2 - 1)2 = 0,0169622, 



and 



J drr 5 K 2 (r)K 2 (r) = 0,08783. 



Putting all this together we have 

«■«. x > = 12^(1 + s> + *$*r ig 2 ' 3 + 5 4 {27,f °' 04874 

+- (2tt) 3 (8 x 0.0169622 - - 0.08783)]. 
4 8 

Working along the same lines with B(a, A) and C(a, A) we find 

, -32A 2 f, f, rp(l-p 2 )Kl{r).Upr) -32A 2 

B(a, A) = . — . / / dp — , ! °\ A J ^ ' = — r^-r x 0,0501125 

v ' ; a 2 (2vr) 3 4! J J F (p 2 + 1) 4 a 2 (2vr) 3 4! 

and 

C < a ' A > = / * ^ = ^1 >< 0. "754499. 

All these numerical values were confirmed by using Fortran. 

Finally, inserting these results in the left hand side of (^|) we get 

1 A A 2 /3 



F(a, A) = — + -— ^ + ——„ - ( 7 + 5.4 x lO" 6 
v ' I2na 2 487r 2 a 2 4\(2ir) 3 a 2 V4 

Comparing this expression with the expansion in A of the right hand side of (^) one 
sees that they are equal up to first order and differ in a small quantity (O(10 -6 )) up to 
second order in A. 

In summary, we have verified the validity of a quantitative prediction of CFT in 
the context of the SG model. Taking into account the well-known bosonization identity 
between the SG theory and the massive Thirring model (characterized by a coupling 
g 2 ) |§ which takes place for /3 2 /(47r) = (1 + g 2 /it)~ l , it becomes apparent that our 
result implies that (|l|) also holds in the strong coupling limit of the Thirring model. 
This, in turn, allows us to improve the proof of the validity of ([!]) for the Baxter and 
Ashkin- Teller models which, up to now, was restricted to the weak coupling limit 
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